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The task of reception of a copy of an arbitrary quantum state with use of 
a minimum quantity of quantum operations is considered. 
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g 1 Introduction 

O I It is known, that an arbitrary quantum state cannot be copied perfectly [1]. However 

^ ' V.Buzek and M.Hillery in work [2] offered a universal quantum cloning machine 

(UQCM), allowing to create 2 identical qubits from 1 qubit. Two output qubits are 



^ . a copy of each other, but they are not a copy of an initial quantum state, and are 

! similar to it only in 5/6 = 0.83. Universality of the UQCM is that it clones any 

qh| quantum state with identical accuracy. In the present work some variants of not 

universal QCM are considered. Their lack is that they clone qubits non-uniformly 
depending on their state. In a final section the calculation scheme of UQCM from the 
general principles is submitted which includes known results [3,4] as special cases. 
It is shown, how one can choose optimum by quantity of used quantum operations 
UQCM from given ones. 

1.1 Brief theoretical data 

Any quantum state 

\^Po) = a\0)+(3\l) (1) 

is represented by a point on Bloch sphere and changes by means of rotation operator 
[5]: 

„ . „ _ f cos 9 —ie ^'-^ sin 9 \ 
^{tl,^) - y sin 9 cos 9 j " 

For simplicity we shall work with equatorial qubits, then = 7r/2 and 

R{9) |0) = cos^|0)+sin^|l), 

R{9) |1) = -sin^lO) + cos0|l) . 
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Pauli matrices can be written down as projectors 

i = ao= |0)(0| + |1)(1|, ai = |l)(0| + |0)(l|, 

a2 = z(|l)(0|-|0)(l|), a3 = |0)(0|-|l)(l|. 
Action of Pauli matrixes on an input is given by expression 

1^^) = iV'o) • (2) 

We shall note at once, that scalar product \ 1IJ2) = 0. It means, that for any 
\iPq) with the help of the universal operator NOT = —ia2 it is possible to create an 
orthogonal state \ip2)- For corresponding Q matrixes of density we shell receive 

pT = \i^^)m■ (3) 

Similarity of two quantum states lip) and |x) is determined by overlapping of their 
wave functions [6]: 

F={^P\p^\^), (4) 

where = \x) {x\- 

The basic two-qubit gate is operation of controlled-NOT (CNOT) [5]: 

P12 \x,y) = \x,xQ) y) (5) 

where © is the modulus-2 summation. The first qubit in this expression refers 
to as control. It does not vary at CNOT transformation. The second qubit is a 
controllable or the target-qubit of CNOT operator. If indexes of the operator to 
change places: P12 = -P21, the second qubit becomes control, and the first becomes 
a target qubit. CNOT operator allows to create the entangled state of two qubits. 
Studying and use of the entangled states is one of the basic problems of quantum 
calculations. 

2 Two qubit QCM 

The ideal cloning of the qubit is forbidden. Really, as Wooters and Zurek have 
shown [1], if there is such a linear unitary operator, that 

U \a) \u) = \a) \a) , 

than for orthogonal states (a| /?) = we have 

U\a + P) \u) = \a + p)\a + P) = \a) \a) + \p) \a) + \a) \P) + \P) \P) . 

On the other hand, owing to linearity U: 

U\a + (3) \u) = U \a) |z/) + U \f3) \u) = \a) \a) + . 
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We came to various expressions that proves impossibility of existence of such op- 
eration. In the present work the les comphcated task is considered: whether it is 
possible to receive, even nonideal copy of any quantum state, and what quantum 
operations are necessary for this purpose? 

2.1 One-operational QCM 

The elementary QCM can be constructed from a single CNOT operator (Fig.l). 
Submitting on one input of the QCM an arbitrary quantum state \ipo), and on other 
- a qubit in state |0) we shall receive 

|v]>-) = |^o)|0) = «|00)+/5|10). 




Figure 1: One-operational QCM 

The work of the QCM is reduced to action by CNOT operator on input qubits 

|^««*) = Pi2|^*"> = a|00)+/3|ll), 

and on an output we receive the entangled state (ebit). The density matrix of 
the output ebit is given by expression = (\Ef°"*|. The reduced density 

operators of the output qubits look like 

pT:^ = a' \0) {0\ + (3' \1) {1\ . (6) 

At the output of the QCM it is received two identical states. It is necessary to 
determine, how output states differ from the input qubit Hjo). Comparing © and 
Oni) it is easy to see, that the output density operator is expressed through the input 
density operator as follows: 

Looking at this expression it is possible to make the assumption, that the target 
state on 50% coincides with input state and has 50% of an impurity. However more 
detailed consideration results in other conclusion. Really, scalar product 

(V'olV'3) = «'-/?Vo, 

i.e. states are not orthogonal, and it means, that wave functions and {ip^) are 
overlapped and part of the information concerning \iIjq) is contained in a density 
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matrix pg". Therefore for calculation of coppying accuracy we shall take advantage 
of the formula (@)) 



(7) 



From last expression it is visible, that accuracy of cloning depends on an input 
state of the original, so received QCM is not universal. It cannot cloning any, 
beforehand unknown states with identical accuracy. Averages d?)) on Bloch sphere 
we shall receive: 



27r 



1 

2^ 



Fde 



3 
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2.2 Two-operator QCM 

Now we shall rotate input |0) qubit for increase of cloning fidelity (Fig.2). Rotation 
of the |0) qubit on angle gives 

IV^^) = /?(0)|O) = cos0|O) + sin0|l). 



|¥o) 







Figure 2: Two-operator QCM 

After that, the state |\E'*") = lipoipji) comes to the input channel of the CNOT. At 
the output of the QCM we'll receive |\Ef°"*) = P12 IvEf*") with the reduced density 
matrixes 



= a' |0) (0| + 2a/5cos0sin0(|O) (1| + |1) (0|) + |1) (1| , 



p°"* = cos2 + sin^ 0) |0) (0| 

+ cos0sin0(|O) (1| + |1) (0|) + (a%in V + /^^ cos» |1) (1| 

and average fidelity 

2 TT 2 1 

Fi = - (cos sin + 1) , F2 = — cos sin + - cos^ + tt sin^ (p. 
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Figure 3: Dependence of the average fidelity Fi 2 on rotation angle for two-operator 
QCM 

As well as in previous case, cloning fidelity depends on the initial state of the 
qubit-original. In addition, accuracy depends on a angle of the rotation |0) qubit. 
Here it is possible to specify some cases (Fig. 3): 

1. at = 0° we shall receive two copies with identical average fidelity: Fi^2 = 3/4. 

2. at = 7r/4 the accuracy dispersion of the first qubit is equal to zero: F^—Ff — 
and the states are separable, i.e. 

|*-*> = p,, iV^o) IV^^) = iV^o) \^r) = |*^"> . 

3. at = 7r/2 measurement results of output states are absolutely anticorrelitive: 

F1F2 — F1F2 _ 
^F^-Ff^Fl-Fi 

4. at = 37r/2 average fidelity of output states again coincides, but we have received 
"no" results: F^ = F2 = 1/2. 

3 Tree-qubit QCM 

More complex case of quantum cloning is carried out by Buzek and Hillery's universal 
QCM (UQCM BH) [2]. Universality of this machine is that its cloning accuracy does 



not depend on a state of the qubit-original. We'll search for the reduced density 
matrixes of output as decomposition on projectors of orthogonal wave functions \iPq) 
and IV'2): 

Pit = foPT + f'pT, (8) 

where /q and /| are so far unknown factors. Substituting in (jS)) decomposition pg" 
and from Q we shall receive 

Pot = i^'fo + (3'f2) |0) {0\+a(3 {fl - fl) (|0) (1| + |1) (0|) + {a'fl + (3^ Q |1) (1| . 

(9) 

It is easy to show, that reduced density Q can be received from three qubits mix 



|M/-*> = |<l>o)oi 10)2 + |$i)oi 11)2, (10) 



•^•o) = n - /I |00) + /3/2 |01) + |10) 



where 



l-^i) = Nfo - fi m + ah |01) + ah |10) 

Then, 



pZ = {c^'fo + |0 ) (0| (11) 
+ 2a(3h^n-n (|0) (1| + |1) (0|) + [a'fl + /?Vo) |1) (1| , 

1 1 

out I e2 e2\ in , r) _c2 in , / e2 e2\ in /-i c\\ 

P2 =2 ~ /2 ) PO + 2/2P1 + 2 (/o - /2 ) P2 • (12) 

In order cross components in © and (lllj) coincided, it is necessary to performe the 
condition 

2/2^ /o - /I = fl ~ fi 

or 

/o=g)/2=g- (13) 

From orthogonality {ipo) and \1IJ2) follows, that the factor of decomposition /q and 
is a cloning fidelity which coincides with average fidelity. 

Sometimes the criterion of universality of quantum cloning is determined by an 
opportunity of representation of the output state as 

pOut ^ ^pin ^ 1^1^ 
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where s is scaling factor. It is easy to show, that the latter expression is a direct 
consequence from orthogonal decomposition (jSI): 

= (/o' - /2) r + ffi- 

Hence there follows connection between the dimensional factor s and factors of 
decomposition /q and /|. 

Small increase of fidelity can be obtained, using phase-covariant (PC) UQCM [3] 
in which the output state is reduced from three qubit mix such as: 

|^°"*> = |0)ol$o)i2 + |l)ol'^i)i2. (14) 

where 

|$o) = ax\00) + py{\Ol) + \10)) + az |11) , 
|$i) = (3z |00) + ay (|01) + |10)) + (3x |11) . 

Then, 



= {a^x' + y') + {y' + z')) \0) {0\ + (15) 
+ 2aP {xy + yz) (|0) (1| + |1) (0|) + {a" {y^ + z") + (3^ {x^ + y^)) |1) (1| . 

Comparing (fTH|) and we see, that in order to obtain the reduced output mix such 
as (jHI) it is necessary to perform the following conditions: 

x^ + y^ = fl (16) 
y' + z' = fl (17) 

2{xy + yz) = fl~fl (18) 

/o+/2=l- (19) 

Obtained optimal task with criterion function (fT^ /q — > max and restrictions (fTTjl- 
(fT^ has the solution 

11 11 1 11 



/2 (max) = - + — =, at x = - + — =, 



z = - 



2 78' 2 Vs' ^ 78' 2 

Let's notice, that if in system (|T6|) - (|T9|) we take z = 0, we shall obtain UQCM BH 

It is obvious, that for the work of the UQCM it is necessary to have in addition to 
the qubit-original two more |00)-qubits. It is necessary to find such transformations 
of input system of qubits \iIjq) ® |0) ® |0), which will result to |\[/°"*) or (fTi)) . 
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It is convenient to divide the process of cloning into two stages (Fig.4). At the 
first stage the entangled state is made from two |00)-qubits with the help of the 
rotate and CNOT operators: 

|00)i2 I^P'-^O = Ci |00) + C2 |01) + C3 |10) + C4 111) . (20) 
For UQCM BH it corresponds to the first part (fTUI) : 

|00)i2 - l^^^'^n = - /I |00) + /2 |01) + h |10) . (21) 
For UQCM PC, in accordance with (fT^ we'll receive 

- l^'^'^O = /o |00) + yi_^(|10) + |01)) + h 111) . (22) 

It is a stage of preparation UQCM for work. At the second stage entangling of input 
iV'o) and prepared qubits by the CNOT operators is made: \iPq) ® |\1>p^'^p) — 

|^f°"*). We shall examine more attentively each of the stages. 



3.1 Preparation of the UQCM for work 

Process of preparation of the UQCM for work consists in obtaining of two |00)-qubits 
of the entangled state with the help of rotation operators R {6) and CNOT: 



IVI/P-P) = i?i(^^3)P2li?2(^^2)Pl2i?l(^^l)|00)i2 (23) 

= (cos 61 COS 62 COS 6'3 + sin 61 sin 62 sin 6^) |00) 

+ (sin 61 cos 62 cos 0-^ — cos 61 sin 62 sin 6^3) |01) 

+ (cos 61 cos 62 sin 6*3 — sin 6*1 sin 62 cos 6^3) 1 10) 

+ (cos 61 sin 6*2 cos 6*3 + sin 61 cos 6^2 sin 6*3) 1 11) . 

Comparing (|20p and (j23|) we'll receive the system 

cos 9i cos 62 cos ^3 + sin 9i sin 62 sin ^3 = 6*1, 

sin 61 cos ^^2 cos ^^3 — cos Oi sin 6*2 sin 6*3 = C2, , x 

cos ^1 cos ^2 sin 6*3 — sin 6*1 sin 6*2 cos 6*3 = C3, 

cos 6*1 sin 0200393 + sin 6*1 cos 6*2 sin ^^3 = C4, 

which has the solution 

/^2 /^2 1 9/^2 9/^2 

2/1 ^2 ^3 2/1 ^^-^2 ^"-^4 

= 1-2CI-2Q + "' 1-201-201 ' 

1 - 2cosW-, ' 
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Now, substituting as the right part values (jUj), we shall have: 
for the UQCM BH 



2. 1/^ 1 



cos^ Qx = cos^ 03 = ^ 

COS 60 = — , 

2 3' 



for the UQCM PC 



cos^^i = cos2^3 = ^l + ^), (25) 
cos^^s = 1. 

Thus, at the first stage we have received angles of rotation operators 

Ri {0i),R2{02),Ri (^3) 

which together with action of CNOT operators allow to create state \'$p^'^p) (pTj) . 
(j22p from two |00)-qubits. We have prepared UQCM to work. Actually, the cloning 
process is carried out at the second stage. 

3.2 A stage of cloning 

Let's choose such combination of CNOT operators which will make an output state 
1^°"*) from input {tpo) and prepared qubits. We'll consider more in detail the 

construction sequence of such operators on the example of UQCM PC. From state 
\ipo) and 



\^P,)0\^^/P^^P) = a 1000) + ^^^^(1001) + 1010)) + /II011) 

+ P (^f^ 1100) + ^t_l (1101) + 1110)) + /2 1111) ) (26) 
we are to obtain the output state (fT^ : 

|^°"*> = «(^/o^|000) + ^^^(|101) + |110)) + /||011) 

+ P (fl |100) + (1001) + 1010)) + fl |111) ) . (27) 



Let's present a state as 

l^o) ® \^r>rep^ = C,y^ |x, y, z) = 1^*"),^, (28) 



and output state p7jl as 



where x, |/, z,p,q,r = 0, 1. 

Comparing values of wave function at identical coefficients Cj, we shall reduce 
transition from to (^Tj) in the table (see Tab.l.) 

Table 1. Truth table for Boolean functions p {x, y,z) ,q (x, y,z) ,r (x, y, z). 



X y z p q r 
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1 





1 
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1 
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1 
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1 








1 








1 





1 
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1 








1 





1 


1 


1 


1 


1 


1 



We have received the table of the validity for Boolean functions of three variables 
p (x, y,z) ,q (x, y,z) ,r (x, y, z). According to the given table we shall write down set 
of disjunctive normal forms and present them as Zhegalkin polynoms: 

p = xSzySzz V xSzySzz V xSzySzz V xSzySzz = x Q) y ® z, 

q = xSzySzz V xSzySzz V xSzySzz V xSzySzz = y, (30) 

r = xSzySzz V xSzySzz V xSzySzz V xSzySzz = z. 

Using the definition of CNOT operator from (j30p we'll obtain 

|*°"*>p,, = \P, q,r) = \x(By(Bz, y, z) = P^^P^^ |x, y, z) = P2A \ v&*">^,,^^ . 

Let's sum up construction UQCM PC (Fig. 4). The first stage is described by the 
expression (f^ . However, according to R2 {62) = 1, therefore (P^j) can be copied 

as 

where the rotation angle of the operator Ri is determined in (I25|) . At the second 
stage two CNOT operators entangle the qubit-original {ipo) with prepared state 
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|o) ] ^ 



prep 



P,7 P7 



Cloning 



Figure 4: Phase-covariant UQCM. After the preparation procedure it is enough to 
work on the qubit-original by two CNOT operators in order to obtain 2 copies with 
accuracy 0.854. 



As result we shall receive an output three qubits density matrix po^g = |\E'°"*) (\['°"*| 
from which with accuracy /o = ^ + — 0.854 two identical bunches = iV'o) (V'ol 
are reduced. 

The similar calculation algoritm at comparison of expressions (fTUI) . (PT|) result 
in sequence of actions of the UQCM BH [7]: 

We shall note the basic differences of two considered UQCM. For UQCM BH the 
pure input state lipo) at the output is registered as a mix (jH)). For UQCM PC the 
input state \iIjq) after cloning procedure has the view 

pT' = IpT + \pT 

and has already 0.25 parts of the impurity. 

The further generalization of calculation algorithms of the UQCM is connected to 
the ability of CNOT operator to interchange the position coefficients Cxyz in decom- 
position ((201). Then, in order to obtain an output state |\['°"*) it is necessary to keep 
only a function |vI/p''^p^ fl20p . and it becomes indifferent where is a factor Cxyz- As an 
example we shall examine UQCM PC ()22|1 . Calculation result are given in Table 2. 
The first line of the second column determines the initial combination of coefficients 
m- (Ci, C2, Cs, C4) = + -ig, ^, -ig, i - -ig) . As C2 = C3, rearranging coeffi- 
cients Cj we shall receive 12 various combinations. For each combination on system 
(1211) ^6 shall find rotation angles 61,62,0^, and on algorithm (j2HI)-(ISni) - operators 
of transition P: |?/'o) ® |\E'P^'^p) — > |\E'°"*). For simplification of record we shall enter 
the following approach: 

arccosJ- + — = - = 22°30', arccosJ- + = = 17°40' 

\/2 78 8 ' y 2 76 540 
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arccos^l^±^ = ^ = 15°, arccos = ^ = 27^20' 

2 12 Y2V 270 

For the description of CNOT operator with inversion the designation is accepted 

-Pi2 y) = Pi2 \x, y) = |x, x e y) = Pi2 \x, aiy) . 

Prom table 2 it is visible, that for each \ there are two ways to obtain |^°"*). It 

is connected to symmetry of an output state with respect to rearrangement {ipi"^^) ^ 
\il^2^^). If we assump that for some combinations (Ci, C2, C3, C4) there are some 
decisions 6*1, 6*2, 6*3 it is quite possible to choose the most simple UQCM from the 
number of obtained UQCM (Fig.4). 
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Table 2. Phase- covariant universal quantum cloning machines. 



N 


(Ci, C2, C3, C4) 


ei 


02 






1 / xyz 




1 


(Ci, C2, C2, C4) 


22°30' 


0° 


22°30' 


xi 


dy®z,y,z) 
dy®z,z,y) 


P12P21P12P10P20 


2 


(Ci, C2, C4, C2) 


27°20' 


15° 


17°40' 


X ( 
X ( 


dz,y,y®z) 
^z,y®z,y) 


P12P20 
P12P21P20 


3 


(Ci, C4, C2, C2) 


17°40' 


15° 


27°20' 


X ( 


^y,z,y®z) 
dy,y®z,z) 


P21P12P10 
P21P10 


4 


(C2, Ci, C2, C4) 


62°40' 


-15° 


17°40' 


X ( 
X ( 


^z,y,y®z) 
dz,y®z,y) 


P12P20 
P12P20P21 


5 


(C2, Ci, C4, C2) 


67°30' 


0° 


22°30' 


x ( 
X ( 


i)y®z,y,z) 
dy®z,z,y) 


P10P20 
P21P12P10P21P20 


6 


(C2, C2, Ci, C4) 


17°40' 


-15° 


62°40' 




dy,z,y®z) 
dy,y®z,z) 


P21P12P10 
PiqPw 


7 


(C2, C2, C4, Ci) 


_17°40' 


75° 


-27°20' 


X ( 


^y,z,y®z) 
dy,y®z,z) 


P21P12PT0 
P21P10 


8 


(C2, C4, C*!, (^2) 


22°30' 


0° 


67°30' 


X ( 

x( 


^y®z,z,y) 
dy®z,y,z) 


P12P21P12PT0P20 
PtqP2Q 


9 


(C2, C4, C2, Ci) 


-27°20' 


75° 


-17°40' 


X ( 

x( 


^z,y®z,y) 
dz,y,y®z) 


P12P20P21 
P12P20 


10 


(C4, C2, C2, Ci) 


67°30' 


0° 


67°30' 


X ( 

X C 


^y®z,z,y) 
dy®z,y,z) 


P12P21 P12 PjoPlo 
P- P- 

-'10-' 20 


11 


(C4, C2, Ci, C2) 


27°20' 


-15° 


72°20' 


x( 
x( 


dz,y®z,y) 
9z,y,y®z) 


P12P21P2O 
PT2P2O 


12 


(C4, Ci, C2, C2) 


72°20' 


-15° 


27°20' 


X ( 
X ( 


^y,y®z,z) 
^y,z,y®z) 


P21P1O 
P21P12P1O 
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